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Abstract 

We study infrared contributions to the QCD description of the HERA data on the structure 
function F2 within the generalized DAS approximation. We argue that this approximation is 
a natural one and consistent with the phenomenon of dynamical mass generation in QCD. The 
investigation is performed at next-to- leading order by using the leading-twist expansion of F2{x, Q^) 
and by adopting an effective charge whose finite infrared behavior is constrained by a dynamical 
gluon mass. We propose one ansatz for the behavior of this effective coupling beyond leading order. 
The dependence of the experimental data on the infrared value of the effective charge is used in 
order to study the asymptotic behavior of the running gluon mass. The deep inelastic structure 
function F2 obtained in this approach shows very good agreement with the experimental data. 
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I. INTRODUCTION 



The nucleoli structure function F2{x, Q"^) at low has been measured in the previously 
unexplored small- a; regime at the HERA collider. The deep-inelastic scattering of leptons 
off nucleons is the instrumental tool for high precision measurements of the quark and gluon 
content of the nucleons and the low transition region bring us into a kinematical region 
where non-perturbative QCD effects becomes essential in order to understand the proton 
constitution. Despite the partonic splitting to be quite well understood through the use of 
the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations and these 
equations being known to describe the data even at not so large Q^, there is no reason to 
expect that they are reliable in the very small-x region. However, perturbative QCD effects 
are expected to become apparent at small-x, where gluon emission off the incoming parton 



leads to power series in as In 



Lipatov (BFKL) equation {2], besides producing an behavior for the gluon distribution, 
generates its own characteristic dependence. Hence approaching the low transition 
region from the perturbative side makes evident the problem of how to incorporate in an 
effective way non-perturbative corrections into the evolution scenario. 

Fortunately, this problem can be properly addressed by bringing up information about 
the infrared properties of QCD, more specifically, by considering the possibility that the non- 
perturbative dynamics of QCD generate an effective gluon mass at very slow region. This 
dynamical gluon mass is intrinsically related to an infrared finite strong coupling constant 
3|, and its existence is strongly supported by recent QCD lattice simulations as well 
as by phenomeno logical results js-?]. It is opportune to remember that phenomenological 
infrared modifications of the strong-coupling constant are quite usual in the literature |8|, |9|, 
nevertheless the fact that an infrared finite coupling constant appears as a consequence of a 
dynamically generated gluon mass is much less known. Furthermore, the dynamical gluon 
mass that constrains the finite coupling constant turns up as the natural infrared cutoff in 
many perturbative QCD calculations, besides being responsible for a smooth transition from 
the perturbative to the non-perturbative QCD behavior 6|, |7|. 

Our task of calculating infrared contributions to the QCD description of the HERA data 
on the deep-inelastic structure function, F2{x,Q^), can succeed in a consistent way by an- 
alyzing exclusively the small-x region since, in this limit, some of the existing analytical 



l/x). Resumming of this series via the Balitsky-Fadin-Kuraev- 



solutions of the DGLAP equation can be directly used 1CH13|. Within this approach the 
HERA data at small- a; is interpreted in terms of the double-asymptotic-scaling (DAS) phe- 



nomenon 



lOj related to the asymptotic behavior of the DGLAP equation in asymptotically 



free field theories 



14l |. The analytical solutions, valid in principle at very small-a; and large- 

jasymptotic part of the evolution, in 



values, can be extended in order to include t 



le su 



17] , leading to small-x asymptotic pre- 



what is called generalized DAS approximation 9|, llSl - 
dictions for parton distribution functions evolved from flat x distributions at some starting 
point Q\ for the DGLAP evolution. In particular, a recent analysis of and its derivatives 
(9^2/(9 In and dXuF^j dXwiX j x) within this approach shows a good agreement with HERA 
data of deep- inelastic scattering for > 1.5 GeV^ 91]. 

One may wonder why the generalized DAS approximation works so beautifully in the 
small-x limit. This fact may be a signal that the choice of a flat distribution for the parton 
distribution function /a(x, Q^) at some initial value Q\ is quite appropriate for the DGLAP 
dynamics, more specifically jisl . 



/a(a;,Qo)=^a {a = q,g) 



(1) 



where Aa are unknown constants to be determined from the data, or in other way, that 
QCD predicts at small-x that F2{x, Q"^) should exhibit double scaling even at not so large-Q^ 
values, provided only that the small-x behavior of the partonic distributions at some initial 
input Qo is sufficiently soft. Actually it was also pointed out that a flat gluon distribution in 
the small-Q^ region appears naturally in QCD [7], within a model for hadronic cross section 
including the phenomenon of dynamical gluon mass generation, which naturally leads to a 
"frozen" infrared effective charge. This mechanism, based on first principles, is probably 
what is behind the good agreement of the generalized DAS approach with the experimental 
data. 

Hence the purpose of this Letter is to compute the structure function F2(x, Q^) of the 



171], assuming the fiat ini- 



proton by means of the generalized DAS approximation 
tial parton distributions as a natural condition for QCD with dynamically generated gluon 
masses, and compare the results with the experimental data of -^2(2;, Q"^) in the infrared 
region. In our calculations the non-perturbative dynamics of QCD is introduced by using 
the infrared finite QCD effective charge dependent on the dynamical gluon mass. As this 
effective strong-coupling has not been determined up to the next-to-leading order (NLO) 



approximation, we propose one ansatz for its behavior at higher order. 

The Letter is organized as follows: in the next section we introduce the generalized 
DAS approach beyond the leading order (LO), and discuss the underlying QCD dynamics 
behind the flat distribution and the frozen effective charge behaviors. In the Sec. Ill we 
propose an ansatz for the NLO behavior of the dynamical strong-coupling based on the 
property of multiplicative renormalizability, showing that this effective coupling reproduces 
the canonical NLO perturbative behavior at large Q^. Our results are presented in the Sec. 
IV, where the analysis of F2{x, Q"^) data is carried on using the formalism developed in the 
previous sections. In Sec. V we present our conclusions. 



II. THE GENERALIZED DAS APPROXIMATION 

The present data of F2(x, Q^) imply a steep gluon at small-x, and there are some successful 
descriptions of F2 by means of DGLAP evolution in the NLO approximation [19]. This steep 
behavior can be generated from a fiat-x gluon distribution at some initial low Qq scale, or 
alternatively it can be directly included into the input distribution to be evolved from some 
higher scale. At sufficiently small- a; we must resum the power series in asln(l/x) via BFKL 
equation. The result of this procedure is sensitive to the infrared fcr region and, for running 
it is found that 

~g{x,kl)^C{kl.)x~\ (2) 

where A ~ 0.5 [20]. Here g{x,k^) is the unintegmted gluon distribution and hence the 
resummation program requires knowledge of the gluon for all k"^ including the infrared 
region. However, in this confinement region the BFKL equation is not expected to be valid. 
Ultimately, with decreasing x, the singular behavior must be suppressed by non-perturbative 
effects. 

The problem of calculating these infrared effects can be addressed by the so called QCD- 
based eikonal models [gI, 2i|, which incorporate soft and semihard processes in the treatment 
of high energy hadron-hadron interactions. At high energies semihard processes are expected 



to give an increasing and significant part of the total hadronic cross sections 22|. Owing 
to the rapid growing of the number of semihard gluons in the hadron at fixed transverse 
momentum, the asymptotic behavior of these cross sections is determined chiefly by the gluon 



distribution. In some QCD-based models these small- a: semihard gluons play a central role, 
and a phenomenological "BFKL- inspired" gluon distribution is introduced [g, l7|, 

gix,Q') = hiQ')x-', (3) 

which captures all the non-perturbative dynamics via the function h{Q'^). Note that here 
g{x,Q'^) is the traditional gluon distribution determined by the parton analysis of the 
F2{x,Q'^) data and whose evolution is controlled by the DGLAP equations, where 



g{x, k 



d{xg{x,Q^)) 



T) 



ainQ2 • 

Recently it was shown that, in order to generate radiatively gluons at small-x, a rapid 
increase of h{Q^) with the momentum has to be accompanied by a fast increase of the J in 
such way that soft values of J are preferred at low [7]. This picture is consistent with the 
statement that the steeply-rising gluon component is absent at low Q"^ and as increases 
it is generated radiatively through perturbative evolution. In the Regge-exchange language 
the quantity J, that controls the asymptotic behavior of the total cross sections, is the 
universal "soft" Pomeron intercept, whose value has been phenomenologically determined 
Q;p(0) = 1 + e ~ 1.1 23|, U^. This is to be contrasted with the "hard" or 

1 + A ~ 1.5. Therefore, fits to a set of hadronic data 



to be J = ap(0) = 1 + e ~ 1.1 
"Lipatov" Pomeron intercept aiiO) 



through QCD-based models show that J starts at a value where the gluon distribution is 
almost flat, J{Q^ ~ ctp(0)5 increases the valence-like character of the gluon 

rapidly disappears 0,0]. It is worth noting that these results are corroborated by a MRST 
analysis of parton distributions of proton 25|. From fitting the sea quark (S) and gluon 



as 



(G) distributions of the defauh MRST partons to the forms fi{x,Q'^) = A{Q)x~^'^'^^^ 
X ^ 0, i = S,G, they have observed that as increases from the input scale Qq = 1 GeV^ 
the fiat behavior of the gluon rapidly disappears due to evolution being driven by the much 
steeper sea. For higher values of the gluon exponent Xq increases rapidly and becomes 
higher than the sea quark exponent A^, since the gluon drives the sea quark via the g qq 
transition. More specifically, A^ starts at a value A^ ~ at ~ 1 GeV^, and by ~ 4 
GeV^ it has the value Ag = 0.2. Hence a fiat input gluon distribution at low momenta, 
that appears as the natural condition for QCD evolution with dynamically generated gluon 
masses, also comes out in standard perturbative procedures. 

In the generalized DAS approach the subasymptotic corrections are included 

via the finite parts of anomalous dimensions of Wilson operators and Wilson coefficients. 



Remember that in the standard DAS approximation only the singular parts of the anoma- 
lous dimensions are taken into account. In the generalized approach the flat input gluon 
distribution ([T]) determines the small-x asymptotics and, at NLO, the twist-two (leading) 
term of F2{x, Q"^) is given by 15l-ll7| 



(5) 



3 47r 

where e = '^{e^/nj is the average charge squared of the effective number n/ of quarks, 
Tr = 1/2 is the color factor for g ^ qq splitting, and 



Ux, Q') = Q') + /-(x, Q') (a = g, g) 



(6) 



the and "— " representation above follows from the solution, at leading twist approxi 



mation, of the DGLAP equation in the Mellin moment space 



16|: 



f^ix, Q') = A-iQ', Qi) exp - + Oix) , 



(7) 



f+{x, Q') = A+{Q', Ql) io{a) exp [-d+{l)s - + 0{p) 



(8) 



/+(x,Q2) = A+(Q^Q,2^ 



1 - — ; 1 phicr) H /o((t) 



47r 



477 



X 



exp[-J+(l)s-D+(l)p] +(9(p) , (9) 
where s = In [as{Ql)/a,{Q% p = [«.(gg) - a,{Q^)] /Att, D±{n) = d±±{n) - {f3Jf3o)d±{n), 



a = 2^ l^d^s + D^pj \nx and P = \j [d+s + Dj^pj /Inx = cr/2 ln(l/a;); here /3o is the 
first (second) coefficient of the QCD /3 function, I^, [v = 0, 1) are functions related to the 
modified Bessel function I^,, and the components of the anomalous dimension d_{n) as well 
as of the singular (d) and regular (d) parts of d^{n) = d^/{n — 1) -\- d^{n), for n — )■ 1, are 



d-i 



AC, 



Ml) 



1 + 



An, 



1-^ 
Cf 



d_(l) 



ACpTLf 



(10) 



/3o 3/3o V ^fJ 3Cyi/3o 

finally, the factors '~ and the components of the singular and regular parts of the remain 
ing anomalous dimensions D± are given by 



A-K 



+ 



C, 



1 - fr 



as{Q\ 
A-K 



diJl 



An 



A 



1 1 



(11) 



A;{Q',Ql)^A,-A^{Q',Ql), 



(12) 



(13) 



A-{Q', Ql) = A, - ^ ^ A^iQ^ Ql) , 



An f 

i++ = -/ (23C^ - 26Cf) , 



20Ca 



(14) 
(15) 



^(36C(3) + 33C(2) 



1643\ 



12 



43 



547 



- 2CK(2) + -Ca-— Cf + 



9 



— — (l-2—\n^ 



36 



-J 

3C| 

2Ca 



(16) 



/So 



4Ci. / , 47 



(17) 



(18) 



20n/ CiT 



di+(l) = 0, 



C.,Ul-2^Jn, 



(19) 



9 Ca' 

where Ca — N , Cp — {N'^ — 1)/2N, and ( is the Riemann zeta function. Prom now on we 
set A?" = 3 in order to fix the Casimir color-factors Ca{— 3) and Cf{— 4/3). With all these 
definitions we can discuss in the next section the QCD effective charge that we shall use in 
our calculation. 



III. THE QCD EFFECTIVE CHARGE AT NLO 

Although not extensively known in phenomenological studies, there is increasing evidence 
that QCD develops an effective, momentum-dependent mass for the gluons, while preserving 
the local SU{3)c invariance of the theory. Of course this mass is not a bare one and at 

7 



few GeV its signal is already erased from the physical amplitudes, which merge into the 
perturbative QCD calculations. In this scenario there is a natural scale that, in principle, 
introduces a threshold for gluons to pop up from the vacuum [26 1. 

Since the gluon mass generation is a purely dynamical effect, the formal tool for tackling 
this non-perturbative problem, in the continuum, is provided by the Schwinger- Dyson equa- 



tions 



271]. These equations constitute an infinite set of coupled non-linear integral equations 



governing the dynamics of all QCD Green's functions. In particular, within this frame- 
work the generation of a dynamical gluon mass is associated with the existence of infrared 
finite solutions for the gluon propagator /S.^y{(f) |28l-l30[|. In covariant gauges, the gluon 
propagator has the form 



(20) 



where ^ is the gauge-fixing parameter and P^uiq) = fi'^i/ — l^i'iu/'f- Infrared finite solutions 
(A~-'^(0) > 0) can be fit by massive propagators on the form A~^(g^) = + m^(g^), where 
m^(g^), which depends non-trivially on the momentum transfer g^, is the so called dynamical 
gluon mass. If the renormalization-group logarithms are included in the Schwinger-Dyson 
analysis, the non-perturbative generalization of the QCD running coupling, the effective 
charge as(g^), is obtained |28l. |29[|. 

Recent studies of a non-linear Schwinger-Dyson equation for t he g luon self-energy show 
that m?{Q'^) may in fact have two distinct asymptotic behaviors 3l|] (note that from now 
on we adopt the virtuality Q in our calculations); first, the dynamical gluon mass runs as an 
inverse power of a logarithm; second, m?{Q'^) drops as an inverse power of momentum. The 
logarithmic running of m'^{Q'^) has been found in previous studies of linearized Schwinger- 



Dyson equations to behave as rn?'{Q'^) ~ (InQ^) 
case 311 this behavior is rewritten as 



-1-7 



with 7 > 1 28|, |32| . In the non- linear 



In 



In 



A2 



-1-71 



(21) 



where 71 = — 6(1 + C2 — ci)/5; here ci and C2 are parameters of the ansatz for the (fully 
dressed) three-gluon vertex used in the numerical analysis of the gluon self-energy. Their 
values are restricted by a "mass condition" which controls the behavior of the dynamical 
mass in the ultraviolet region. In the case of a logarithmic running, ci G [0.15,0.4] and 
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C2 G [—1.07, —0.92] [31]; the values of the parameters p and m^, which control the behavior 
of m^lQ"^) in the infrared region, are also restricted by the mass condition, and general con- 
straints are satisfied for p G [1.0,8.0] and nig G [300,800] MeV |33|. It is worth mentioning 
that the dynamical gluon mass was found for the first time by Cornwall to be equal to 



ln(^ 



12 
' 11 



(22) 



where the infrared mass value is phenomenologically determined and typically of the 



order m„ = 500 ±200 MeV h 



m 



28 



34| . Note that the Cornwall expression f l2^ is a special 
case of the logarithmic running mass ( 12 ip : more specifically, f l22p can be obtained from ( 12 ip 
by fixing p = 4 and 71 = 1/11. 



A power-law running behavior for iv?{Q'^) was first envisaged in 28|, |35|]. According to 
an OPE calculation the most probable asymptotic behavior of the running gluon mass is 
proportional to 1/Q^ 36|. At the level of an non- linear Schwinger-Dyson equation this 
asymptotic behavior is given by 



m2(Q2 



ml 



mi 



In 



A2 



In 



is? 



(23) 



where 72 = (4 + 6ci)/5; for power law running the mass condition imposes ci G [0.7, 1.3]; 
the p and parameters are constrained to lie in the same interval as before, namely 
p G [1.0,8.0] and mg gJSOO, 800] MeV 



31 



The results of Ref. 



33] 



3l| are precise with respect to the gross asymptotic behavior of the 



running gluon mass which are represented by our equations (21) and (23). Eq.(5.2) of that 
reference contains a broad definition of the mass function used to match the numerical re- 
sults, and it must be said that the approximations in there exclude the ghost fields and 
the regularization of quadratic divergences is obtained through what is called "tadpole- 
condition". This procedure does not determine a;s(0) and mg univocally, and leads to a 
dispersion on the frozen coupling infrared behavior. However, this uncertainty is system- 
atically reduced by the aforementioned phenomenological studies ^-7], which determine a 
frozen value for the LO effective charge of the order q;s(0) ~ 0.7 ± 0.2. Our result for the 
NLO frozen behavior, a^^^{^) ~ 0.6, despite not being directly comparable to LO results, 
gives support to the statement that the dynamical gluon mass mg is not strongly dependent 



on the perturbation order. Moreover, our value 0^ (0) ~ 0.6 is totally consistent with the 



frozen value as(0)/7r ~ 0.19, obtained very recently from an analytic QCD model 37l |. 

Given the running behavior of m'^{Q'^), the leading-order QCD effective charge asiQ"^) is 
written as 

aJQ^) = -J (24) 



where Bq = /3o/47r = (l/47r)[(llCA — 2nj)/3] and A = Aq^^. The effective charge clearly 
shows the existence of an infrared fixed point as 0, i.e., the dynamical mass term tames 

the Landau pole and as freezes at a finite value in the infrared limit. It must be stressed 
that the fixed point does not depend on a specific process, it is uniquely obtained as we fix 
A and, in principle, it should be exactly determined if we knew how to solve QCD. Note 
that in the limit ^ A^ the dynamical mass m{Q'^) vanishes, and the effective charge ( 12^ 
matches with the one- loop perturbative QCD coupling ^^(Q^). It means that the asymptotic 
ultraviolet behavior of the LO running coupling, obtained from the renormalization group 
equation perturbation theory, 

af^(Q^ » A^) ^ , (25) 

is reproduced in solutions of Schwinger-Dyson equations, provided only that the truncation 
method employed in the analysis preserves the multiplicative renormalizability (MR). Since 
the MR is an important feature of gauge field theories, and holds for any renormalization 
scale, we argue that a QCD effective charge at NLO, a^^^ , can be successfully built by 
saturating the two- loop perturbative strong coupling a^^'-^, that is, by introducing the 
replacement af^'^(Q^) af^^iQ"^) = af^*^(Q^ + Am^iQ^)) into the perturbative result. 
In this way, the QCD effective charge at NLO is given by 

In (in (Q!±^!!l!m^ 



6oln(^^^±^ 



(26) 



where hi = Pi/lGir^ = (l/167r2)[(34Cl - ^^(lOC^ + 6C^))/3] and A = A^^g. Note that in 
the limit ^ A^ the effective charge (l26ll matches with the canonical two-loop perturbative 
coupling, a^^*^, in such a way that the relation 



^NLO ^NLO 



(27) 
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is valid in the ultraviolet region. This relation is expected to be valid if the Schwinger- Dyson 
equation is renormalized multiplicatively. 

We have created other ansatzes preserving the relation (!27l) . where an intermediate scale 
was introduced in order to separate the perturbative and non-perturbative regions, but they 
did not introduce significant differences in the behavior of ^« in the infrared region. Thus 
we have adopted the coupling ( 126|) as the standard ansatz in our calculations of the structure 
function F2{x,Q'^) of the proton. We do not expect that any phenomenological calculation 
using the logarithmic and power-law running dynamical masses will be strongly dependent 
on the asymptotic behavior as they are on the infrared one. Therefore, since the calculation 
is quite dependent on the behavior of the effective charge in the infrared region, and our 
analysis includes HERA data sets at low and moderate Q^, we carry out two independent 
global fits to HERA data: in the first one we adopt the effective charge (l26l) with a loga- 
rithmic mass running, expression (!2T|) : in the second fit the effective charge runs through 
the power-law mass running, expression (l23l) . Within this procedure we can investigate if 
the experimental data can differentiate these solutions. The different momentum behaviors 
of the canonical (perturbative) and the effective charges are shown in FigJH 



IV. RESULTS 



From the formalism discussed in the previous sections, we analyze F2{x, Q ) data sets at 



low and moderate values 38| , by adding the statistic and systematic errors in quadrature. 
We carry out global fits to F2 data by means of a fitting procedure with an interval ~ 
Xmin corresponding to the projection of the hypersurface containing 90% of probability. 
To keep the analysis as simple as possible, we fix n/ = 4 and A = 284 GeV. These choices 
are not only consistent to NLO procedures, but are also the same ones adopted in Ref. jo]. 
Concerning the QCD effective charges, we fix in all the fits p = 4, since this is the optimal 
value obtained by Cornwall in order to reproduce the numerical results of a gauge invariant 
Schwinger- Dyson equation for the gluon propagator [28]. Moreover, we observe that the 
gluon mass scale rUg is not very sensitive to p (at least in this leading-twist operator analysis), 
changing by about 12% (4%) when p ranges from 1.0 to 8.0 in the case of a logarithmic 
(power-law) running behavior. 

Our first analysis consisted in the determination of the Ag, Ag and Ql values from a global 

11 



fit to F2 data using the canonical (perturbative) QCD coupling at NLO. The /DoF for this 
fit was 2.88. These values are shown in Tabled and the structure functions corresponding to 
these values are shown by the dotted-dashed curves in Fig|2J It is clear from the relatively 
high value of 'x^ / DoF obtained in this fit, as well as from the curves depicted in Fig|21 that 
the canonical version of the coupling constant provides a worse fit for the x dependence 
of F2{x,Q'^) for specific bins, and the disagreement is larger for smaller values of Q^. 
This result was expected in the light of an earlier analysis on F2 data, which showed the 
requirement for theoretical improvements in the QCD canonical coupling for smaller 
values 9j. More specifically, it was shown that modifications in the strong coupling based on 
the Kallen-Lehmann analyticity 39| , or on a purely phenomenological freezing, improve 
the description of the experimental F2 data at low Q^. Our results for F2{x,Q'^) in the 
case of the standard perturbative coupling are therefore similar to the ones ofRef.M. It 

CI 

is important to stress that the analytic coupling constant used in the Ref.[9| has a frozen 
value as{0) 1.398, which is not consistent with phenomenological results using the gluon 
dynamical mass. Moreover, this analytic coupling, as pointed out by Cvetic, Kogerler and 



Valenzuela 37|], does not give the correct value of the well- measured semihadronic r d ecay 



ratio r^-, namely r'f^^ = 0.203 ±0.004. However, the analytic coupling version of the Ref. 37|, 
besides generating a frozen value similar to the one obtained in this Letter, reproduces 
successfully the experimental value of Vr- 

In the second analysis we carried out a global fit to F2 data using the QCD effective 
charge with a logarithmic running mass, namely, using the expressions fl2Tl) and fl26l) . In 
principle, the infrared mass scale rrig and the factor 71 are also, together with Ag, Ag and 
Qq, fitting parameters to be determined, where the constraints on ci and C2 are satisfied for 
7i e [0.084,0.564]. However, the best value for this fit is obtained at 71 = 0.084. If we 
carry out global fits to different combinations of F2 data sets, as for example taking into 
account only data sets with < 2.0 GeV^ or > 2.0, we observe that the optimum 
for each fit has all 71 = 0.084. Therefore, we have set 71 = 0.084 in our subsequent analysis. 
The x^ / DoF obtained by this fit was 1.87. The values of m^, Ag, Ag and Qq are show in the 
second line of the Table [H The structure functions corresponding to these values are shown 
by the solid curves in FigjSJ 

In the sequence we have obtained rUg, Ag, Aq and Qq by means of a fit to F2 data using 
the QCD effective charge with a power-law running mass, namely, using the expressions 
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(123|) and ( 126|) . In this analysis 72, constrained by C2 to lie in the interval [1.64,2.36], is set 
to 72 = 2.36, since the optimum for fits to different combinations of data is obtained 
for 72 = 2.36. The x^/DoF for this global fit was 2.13. The values of m^, Ag^ Ag and Qq 
are show in the third line of the Table [H The theoretical F2 results corresponding to these 
values are shown by the dashed curves in FigJ21 

Note that all Qq dependence enters into the definition. In our approach Qq will always 
appear added to a factor 4m^((5o)- Particularly in the case of log-running gluon mass 
the Ql dependence is not felt strongly, as a consequence of a very fiat infrared coupling 
constant. In this case the "effective" initial scale Qe//o dominated by the term 4m^(Qo)) 
namely Qe//o ~ 0.009 + 4m^(0. 009) ~ 0.534 GeV^. On the other hand a power-law running 
gluon mass leads to a ag behavior that matches very fast with the standard perturbative 
one and consequently a stronger dependence on the Ql parameter, where Qlff^ = 0.029 -|- 
4m2 (0.029) ^ 0.456 GeV^. These values are close to the result from Ref.Q] for the NLO fit 
using a frozen coupling, namely Ql = 0.589 ± 0.006 GeV^. 

It seems that owing to the larger x^/DoF value obtained in the power-law mass analysis 
the data favors the logarithmic running mass case, but it must be said that the differences 
are quite subtle and in both cases we have a substantial improvement in the agreement of 
the theory with the experimental HERA data at low Q^. 



V. CONCLUSIONS 

In this Letter we have computed the structure function F2{x, Q"^) of the proton by means 
of the generalized DAS approximation with a QCD effective charge at NLO. This effective 
strong coupling is finite in the infrared region and naturally connected to the phenomenon 
of dynamical gluon mass generation in QCD. Its basically fiat behavior below the rrig scale 
indicates the existence of an infrared fixed point. 

We have observed that an infrared finite coupling constant is fundamental in order to 
improve the description of the F2{x, Q"^) experimental data at low Q^, what can easily be seen 
if we compare the results using effective charges with the one using the canonical perturbative 
coupling. Such fact was already observed in Ref.jsjl, in which a phenomenological infrared 
finite coupling was also adopted. We stress that our QCD effective charges, despite the same 
freezing behavior of the couplings adopted in j^j], are obtained from the QCD Lagrangian, 
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i.e., they are derived from first principles, in a scenario wliere infrared effects are taken 
into account. Moreover, we have argued that a flat x initial condition in the DGLAP 
evolution equations, which determine the basic role of the singular parts of the anomalous 
dimensions in the generalized DAS approach, is naturally related to what is expected from 
non-perturbative QCD, i.e., from QCD with dynamically generated gluon masses. 

Through global fits to F2 data we have obtained the best values of the infrared mass 
scale in the case of the logarithmic and power-law running mass, rrig = 364 ± 26 MeV and 
TTig = 355 ± 27 MeV, respectively. It is important noting that these infrared scale values 



are of the same order o: 
interacting processes 



magnitude as the values obtained in other calculations of strongly 
-7|. These results corroborate theoretical analyzes considering the 
generation of a dynamical gluon mass in non-perturbative QCD. 

Our results show that the leading-twist approximation of the Wilson operator product 
expansion is quite accurate on the description of the structure function F2 data. However, 
our results using QCD effective charges indicates that, at principle, in order to differentiate 
the logarithmic running mass from the power-law one, a higher-twist study is necessary. An 
analysis using higher twist corrections to the expansion of F2{x, Q^) is in progress. 
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FIG. 1. The canonical (perturbative) coupling constant and the QCD effective charge with loga- 
rithmic and power-law mass running at NLO. 
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FIG. 2. Fits of the x dependence of F2{x,Q'^) for specific with the QCD effective charge 
dependent on a dynamical gluon mass. The sohd curves were obtained considering the logarithmic 
running gluon mass, Eq. (|2ip . whereas the dashed curves correspond to the power-law running case, 
Eq. (p3|) . The dotted-dashed curves correspond to the perturbative behavior at NLO. 
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Coupling 


nig [MeV] 


A, 


A, 


Ql [GeV2] 


xVDoF 


Canonical 




-0.339±0.019 


1.119±0.025 


0.414±0.016 


2.88 


Logarithmic 


364±26 


-0.084±0.063 


0.843±0.069 


0.009±0.116 


1.87 


Pow(M;-Law 


35o±27 


-0.253±0.011 


1.018±0.5G5 


0.029±0.008 


2.13 



TABLE L Values of the parameters rUg, Ag, Aq and Qq resulting from the global fit to F2 data. 
The errors were obtained assuming a confidence region of the parameters of 90%. 
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